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Abstract -In this paper, the hydrodynamics of a direct contact heat transfer in a Bubble Column Reactor 
are investigated using a three dimensional, Eulerian-Eulerian Computational Fluid Dynamics (CFD) model. In 
particular, the effect of the static liquid (H) height on gas holdup (aş) is examined by using three different heights 
(45, 55, and 65 cm). Using experimental data, the model is validated and compared to a previous two-dimensional 
CFD model. It is found that the three-dimensional system can precisely model the trends of gas holdup as the 
static liquid height height changes. In addition, compared to the 2D model, the 3D model is more precise and 
overestimated the gas holdup. In addition, increasing the static liquid height causes the gas holdup to decrease; 
where the gas holdup decreases by 11% when the static liquid height increases from 45cm to 65cm. 
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1. Introduction 

It is anticipated that hydrogen will be a major contributingcomponent to the sustainable energy 
supply in the future [1] since the usage of hydrogen will minimize the pollution that contributes to 
climate change by lowering greenhousegas emissions. Hydrogen gas generation from a variety of 
fuel gases may create greenhouse gases, however the emissions are much lower than those typical 
produced by gasoline and diesel vehicles. Thermochemical cycles are prospective solutions that can 
be coupled with nuclear reactors in order to thermally decompose water into oxygen and hydrogen 
via multistage processes. The copper-chlorine (Cu-Cl) cycle has been identified by Argonne 
National Laboratories (ANL) as one of the most promising low temperature cycles [2], [3]. 

The oxygen reaction in the Cu-Cl cycle of hydrogen production is a high temperature reaction 
that requires a source of high temperature heat. This heat may be generated using nuclear reactors or 
solar thermal energy which are non-polluting source of high-temperature heat. It has been found that 
the more practical and effective method of heating oxygen reactor is to heat the molten salt inside the 
oxygen reactor, which will transmit heat from the molten CuCl to the solid Cu2OCI2 (reactant) 
particles within the reactor. Different mechanisms of heat transfer for the oxygen reactor have been 
investigated [4-8]. It has been determined that direct contact heat transfer from the oxygen gas to the 
molten CuCl is the optimal heat transfer mechanism for the oxygen reactor [8-9]. In this method, a 
portion of the oxygen gas generated by the oxygen reactor's decomposition process is heated to 530°C 
and reinjected into the oxygen reactor to transfer heat directly to the molten salt. 

Research on bubble column reactors has been done both experimentally and numerically 
throughout the years. The Eulerian approach to CFD analyses for bubble column reactors has been 
reviewed by Abdulrahman et al. [10]. Abdulrahman [11—15] has examined the results of 2D CFD 
simulations for the bubble column reactor's gas holdup, volumetric heat transfer coefficient, gas and 
slurry temperatures, and solid concentration. Sarhan et al. [16] have investigated the effects of the 
physico-chemical properties of the liquid and gas phases on bubble formation and hydrodynamics of a 
bubble column reactor using the population balance equation combined with a 3D CFD model. They 
have used a Euler-Euler CFD model to predict experimental results of the gas holdup in a bubble 
column reactor using different phase flows within the range of +7%. Also, they have observed that the 
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gas holdup will increase slightly as the gas phase density increases. Li et al. [17] have conducted both 
experimental and CFD hydrodynamic analysis of an air-water- glass beads slurry bubble column. The 
slurry bubble column reactor was modeled using a 2D axisymmetric two fluid Euler k-€ model. It was 
noted that the change of hydrodynamic characteristics with column diameters is the major cause of 
bubble column scale-up rules. Bubble columns with wider reactor diameters result in the axialliquid 
velocity rising dramatically within the column core, while the gas holdup is very minimally 
influenced. Also, it was noted that with increasing column scales, turbulent kinetic energy rises. 

Li and Zhong [18] conducted a 3D, Eulerian-Eulerian-Eulerian, three phase (air-water—glass 
powder), time dependent, CFD analysis of three different bubble column reactors to studythe 
hydrodynamics in relation to time step, momentum discretization schemes and wall boundary 
conditions. The three different models used where the Gandhi et al. (height:2500mm, Static height: 
1500mm diameter:150mm),the Rapure et al. model (height: 2000mm static height: 1000mm 
diameter:200mm) and the Li and Zhong model (height:800mm width: 100mm depth:10mm). The 
turbulence model used for their simulations was the RNG k-€ model. Li and Zhong concluded that 
the conditions that best reflected experimental results were the use of a no slip condition, momentum 
discretization using the second order upwind, and a time step of 0.001s. Pu et al. [19] conducted a 
2D CFD simulation of a molten salt bubble column, to investigate the hydrodynamics and direct 
heat transfer characteristics of a two-phase flow model (air —molten salt). An Euler- Euler 
multiphase model with a k-€ turbulence model was used. Factors investigated during the simulation 
include changing the superficial gas velocity, varying the static liquid heights, and using different 
operationalpressure and inlet gas temperature. It was observed that as the superficial gas velocity 
or the operational pressure increases so does the molten salt temperature and rising rate of the 
molten salt temperature over time. However, when the static liquid height rises, the rate of increase 
in average molten salt temperature falls. Increases in superficial gas velocity or operating pressure 
enhance the volumetric heat transfer coefficient, whereas increases in static liquid height lowersthe 
volumetric heat transfer coefficient. 

Zhang and Luo [20] generated a 2 phase (air-water) CFD model to investigate a bubble 
column’s local gas-liquid slip velocity distribution in relation to heat transfer in a heterogeneous 
regime. The study also investigated the simulated time average of the local 2 phase slip velocities 
when varying superficial gas velocities, axial locations and scale of the bubble column. The model 
used was a CFD -PBM (population balance model) simulation with a RNG k-€ turbulence. It was 
observed that raising the superficial gas velocities raised the local gas-liquid slip velocities in the 
region of developed flow. At larger superficial gas velocities, it was notices that the slip velocities 
were more affected by the radial position. The slip velocities near the center of the column where 
lower than that of the fully developed region. The slip velocities were minimally affected by the 
axial heights for the fully developed flow regions. Li et al [21] investigated the effects of a circular 
heat exchanger using a 2D CFDPBM model on the hydrodynamics of a pilot scale slurry bubble 
column reactor. AEuler —Euler multiphase model with an RNG k-€ turbulence model was used. The 
reactor investigated had a diameter of 30cm, height of 200cm and the circular heat exchanger had 
a height of 108cm. Paraffin oil and catalyst particles where the assumed materials in the simulation. 
It was observed that the gas phase was notably distributed, local circular vertices were generated, 
and the slurry was strongly circulated due to the implementation of the circular heat exchanger tube. 
The bimodal profile of the gas holdup profile in the radial direction is caused by the circular gas 
distributor's particular layout. Furthermore, the circular heat exchanger tube increases this 
distribution, resulting in a larger gas holdup, which facilitates momentum transfer. Based on the 
above literature review, it can be concluded that the effect of static liquid height on gas holdup in 
the oxygen bubble column reactor of the Cu-Cl cycle for hydrogen production using 3D CFD 
simulations has not previously been investigated. This paper addresses the aforementioned gap. 


2. CFD Simulation Model 

The simulations of this thesis are validated against the experimental data provided by Abdulrahman 
[9, 22-25]. As such the simulated reactor is designed to be a simplified version of the physical reactor 
to reduce computational costs. The experimental reactor was constructed of stainless steel with a 
diameter of 21.6 cm. A stainless-steel distributor is inserted 10.8cm above the base of the SBCR. The 
gas was fed into the SBCR using a six-arm sparger type gas distributer. Each arm of the sparger had 12 
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orifices with 0.3cm diameters (72 holes in total). The orifices on the sparger were 4.4 cm, 5.5, 6, and 
6.8cm from the center of the reactor [4] as seen in Fig. 1. 


Bottom View Side Sectinal View 


Fig. 1 Sparger design [9]. 


Initially the reactor and sparger designs were recreated using Inventor Professional. However, 
in order to reduce computational costs, the designs were simplified. First the sparger head was 
lowered by 10 cm so that the gas would be released at the base of the reactor. This would mean that 
the simulated heights (H) would be comparable to the experimental results that are 10cm higher. 
For example, the reactor heights in the simulations are 35, 45 and 55cm which were compared to 
the experimental results when the static liquid heights were 45, 55, and 65cm. The model could be 
further simplified to reduce computational costs. As such the sparger was simplified to a single inlet 
at the base of the reactor opposed to 72 inlets. The inlet has a diameter of 18cm and is extruded 
0.3cm downward. The diameter for the gas inlet was selected to represent the bubble distribution at 
the base of the experimental reactor. The progression of the design simplification can be viewed in 
Fig. 2. The final reactor design can be seen in Fig. 3 with overall diameter 21.6 cm, gas inlet 
diameter18cm extruded 0.3cm downward with varying heights depending on the experiment. 
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(a) (b) (c) 
Fig. 2 Progression of simplifying the reactor design (a) most complex, (b) more simplified (c) most simple (final 
design). 


0.3 cm | y 


18 cm 
21.6 cm 
Fig. 3 Final SBCR geometry used for simulations. 


Due to the difficulties associated with the Cuprous Chloride (CuCl) and Oxygen (O2) substances, 
Abdulrahman [9, 22-25] conducted his experiments on a Water-Helium-Alumina system. The 
challenges include the difficulty in observing O2 bubbles in molten CuCl because of its dark colour, 
the corrosiveness of the CuCl molten salt, and the O2 gas's ability to oxidize a variety of substances to 
accelerate its combustion [9, 26-27]. Based on the Buckingham pi theorem, a dimensional analysis 
revealed that liquid water at 22 °C and Helium gas at 90 °C are suitable substitutes for molten CuCl at 
530 °C and oxygen gas at 600 °C [9]. In this paper, simulations of a 3D plane system are performed 
using an Eulerian-Eulerian model, an Eulerian sub-model, and a pressure-based solver type. The 
turbulence model used is RNG k-€ model and the wall function is the standard wall function. For the 
BCR, a hexahedron mesh is used. The independence analyses of the mesh are performed to guarantee 
that the largest mesh size is chosen to minimize computational costs while achieving acceptable 
results. The mesh's final composition consists of 26,825 nodes and 24,396 elements. This resulted in a 
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3% difference in gas holdup when utilizing finer meshes. The inlet boundary condition is specified 
with an inlet superficial gas velocity and a gas holdup of 1. The pressure at the outlet is set to 
atmospheric pressure. For both phases, a no-slip condition is applied to the reactor's walls. It is 
challenging to estimate turbulent kinetic energy and dissipation rate at the inlet and outlet, so 5,000 
iterations were used to specify them. Table 1 displays the equations used in the CFD analysis of this 
paper. The equations in Table 1 are written exclusively for the gas phase. Since the equations for 
liquid phase are similar to those for gas phase, they are not repeated. 


Table 1: Details of equations used in the 3D CFD simulations. 


Description [reference] Equation 
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3. Results 

Figure 4 shows the three-dimensional gas holdup (æg) versus static liquid height (H) and 
superficial gas velocity (Ugs) curves. The effects of varying H on æg while changing the superficial 
gas velocity (Ugs) for a helium-water system are shown in Fig. 5. It can be observed from the graphs 
that increasing H decreases æg. The contours of the cut section of the BCR taken in the center of 
the XY and ZY planes are shown in Figs. 6 and 7. Additional cut sections are taken at varying 
heights on the ZX plane within the reactor at heights of 10, 20, and 30 cm from the reactor's base 
in order to obtain a more precise profile of æg. It is evident from the contours that the gas holdup is 
not symmetrical in the XY, ZY, and ZX planes, indicating that its behaviour is highly three- 
dimensional. It has been observed that as H increases, ag decreases. At a superficial gas velocity of 
0.05m/s, increasing H from 45cm to 65cm reduces the gas holdup by approximately 11%. With a 
superficial gas velocity of 0.15m/s, an increase in H from 45cm to 65cm resulted in a decrease of 
approximately 15% in the gas holdup. This is due to the increase in hydrostatic pressure and the 
decrease in pressure as H increases for specific Ugs. Additionally, having a shorter BCR (Hr/Dr < 
3) can prevent the complete development of liquid circulation, thereby decreasing æg [30]. 
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0.15 ™ 25-30 
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m 15-20 
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oo 
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= 55 0.05 


H (cm) 65 
Fig. 4 Average Gas Holdup Versus Static Liquid Height and Superficial Gas Velocity of Helium-Water system. 
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Fig. 5 The effect of static liquid height on gas holdup for different superficial gas velocities. 
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Fig. 6 Gas holdup contours for different static liquid heights and Ugs = 0.05 m/s. 
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Fig. 7 Gas holdup contours for different static liquid heights and Ugs = 0.15 m/s. 
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The gas holdup results of the 3D-CFD Helium-Water BC simulations were validated against 
experimental data from Abdulrahman [9] (Fig. 8) and compared to 2D simulations from Abdulrahman 
[9] (Fig. 9). Figures 8 and 9 show the ability of the 3D-CFD models to determine the gas holdup (æg) 
at various static liquid heights (H). Theoretical 3D-CFD models can predict the experimental data 
with reasonable accuracy. It can be observed that the majority of gas holdup simulation results were 
overestimated. Reducing the mesh size is one possible technique for minimizing relative error. The 
software will be able to account for minor vortical structures in the flow, such as eddies, if the mesh 
size is decreased [31]. 


35 Exp Ugs=0.05m/s 


Exp Ugs=0.08m/s 
Exp Ugs=0.1m/s 
Exp Ugs=0.13m/s 


e E è p> x 


Exp Ugs=0.15m/s 
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----- 3D CFD Ugs=0.1m/s 
seeeeeeee 3D CFD Ugs=0.13m/s 


3D CFD Ugs=0.15m/s 


45 50 55 60 65 
H (cm) 
Fig. 8 Comparison of the average gas holdup versus static liquid height between 3D CFD simulations and 
experimental data. 
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Fig. 9 Comparisons of the average gas holdup versus superficial gas velocity for different static liquid heights 
between 3D CFD simulations, 2D CFD simulations and experimental data. 


5. Conclusions 

This paper investigates 3D CFD analyses for the hydrodynamics of a bubble column reactor with 
direct contact heat transfer. The results of this paper demonstrates that gas holdup flow patterns in 3D 
CFD simulations are asymmetrical. Also, gas holdup increases as liquid height decreases. This is due 
to the increase in hydrostatic pressure and pressure drop that occurs when the static liquid height of 
the reactor increases at specific superficial gas velocity. With a superficial gas velocity of 0.05m/s, it 
is observed that increasing H from 45cm to 65cm decreases the gas holdup by 11%. Moreover, it is 
observed that the 3D-CFD results slightly overpredict the experimental results for gas holdup, 
whereas 2D simulations underpredict the experimental results. Finally, it is possible to conclude that 
3D simulations are more accurate than 2D-CEFD simulations. 
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List of Symbols 


A; Interfacial area concentration Velocity field 
C Specific heat Vg Volumes of gas 
Cp Drag coefficient VI Volumes of liquid 
dp Bubble diameter Ugs Superficial gas velocity 
g Gravitational acceleration ag Gas holdup 
H Height Here Effective viscosity 
. Total interfacial forces between the Dynamic viscosity gas 
Mi H 
phases g 
P Phase pressure Ly Dynamic viscosity liquid 
Q Intensity of heat exchange between the Pg Density, gas 
2l gas and liquid phases 
Re Reynolds number Pl Density, liquid 
T:VV Viscous stress tensor contracted with 
T Temperature 


the velocity gradient 
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